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ABSTRACT 

The influence of tlie metallicity at tlie main sequence on tfie cfiemical structure of the 
exploding white dwarf, the nucleosynthesis during the explosion and the light curves of 
an individual Type la supernovae have been studied. Detailed calculations of the stellar 
evolution, the explosion, and light curves of delayed detonation models are presented. 

Detailed stellar evolution calculations with a main sequence mass Mms of 7 Mq have 
been performed to test the influence of the metallicity Z on the structure of the progen- 
itor. A change of Z influences the central helium burning and, consequently, the size of 
the C/0 core which becomes a C/0 white dwarf and its C/0 ratio. Subsequently, the 
white dwarf may grow to the Chandrasekhar mass and explode as a Type la supernovae. 
Consequently, the C/0 structure of the exploding white dwarf depends on Z. Since C 
and O are the fuel for the thermonuclear explosion, Z indirectly changes the energetics 
of the explosion. 

In our example, changing Z from Pop 1 to Pop 11 causes a change in the zero point 
of the maximum brightness/decline relation by about 0.1'" and a change in the rise 
time by about 1 day. Combined with previous studies, the offset in the maximum 
brightness/decline AM ?a At where At is the change of the rise time in days. 

Systematic effects of the size dissussed here may well make the results from the 
SNe la searches consistent with an Universe with = 0.2 and I^a = but hardly 
will change the conclusion that we live in a universe with low ^m- Variations of the 
expected size may prove to be critical if, in the future, SNe la are used to measure large 
scale scalar-fields because Z may show large local variations. Evolutionary effects will 
not change substantially the counting rates for SNe la even at very large red-shifts. 

Evolutionary effects may be of the same order as the brightness changes related to 
cosmological parameters, but we have shown ways, how the effects of evolution can be 
detected. 

Subject headings: Supernovae: general - nucleosynthesis - radiation transfer - large-scale 
structure of universe 
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1. Introduction 

Two of the important new developments in obser- 
vational supernova research in the last few years were 
to establish the long-suspected correlation between 
the peak brightness of SNe la and their rate of decline 
by means of modern CCD photometry (Phillips 1993) 
and the exact distance calibrations provided by a HST 
project (e.g. Saha et al. 1997). This allowed an em- 
pirical determination of Ho with unpreceded accuracy 
(Hamuy et al. 1996, Riess et al. 1995). Independent 
from these calibrations and empirical relations, Ho 
has been determined by comparison of detailed theo- 
retical models for light curves and spectra with obser- 
vations (Miiller & H5flich 1994, Hoflich & Khokhlov 
1996, Nugent et al. 1996). All determinations of the 
Hubble constant are in good agreement with one an- 
other. More recently, the routine successful detection 
of supernovae at large redshifts, z (e.g. Perlmutter et 
al. 1999, Riess et al. 1998), has provided an exciting 
new tool to probe cosmology. This work has provided 
results that are consistent with a low matter density 
in the Universe and, most intriguing of all, yielded 
hints for a positive cosmological constant. These re- 
sults are based on empirical brightness-decline rela- 
tions which are calibrated locally. This leaves sys- 
tematic effects as the main source of concern. In this 
respect, theoretical models provide a critical tool to 
estimate the possible size of evolutionary effects, how 
these effects can be recognized and how one may be 
able to correct for them. 

The primary scenario for SNIe la consists of mas- 
sive carbon-oxygen white dwarfs (WDs) in a close bi- 
nary system which accrete through Roche-lobe over- 
flow from low mass companion star when it evolving 
away from the main sequence or during its red giant 
phase (Nomoto & Sugimoto 1977). The WD is the 
final evolution of stars with a main sequence masses 
smaller than « 8 Mq which has lost his H/He rich 
envelope with the C/0 core of « O.5...1.2M0. If ac- 
cretion is sufficient large (i.e. > 2. ..4 lQ~^MQ/yr), 
the accreted H burns to He and, subsequently, to 
C/0 on the surface of the WD and the mass of the 
WD grows close to the Chandrasekhar mass Mch- 
The corresponding time scales for growing to Mch 
are < 5 x IQ^yrs. In these accretion models, the ex- 
plosion is triggered by compressional heating close to 
the center. 

Hoffich, Khokhlov & Wheeler (1995) showed that 
models based on Chandrasekhar mass carbon-oxygen 



white dwarfs can account for subluminous as wells 
as "normally bright" SNe la. The basic paradigm 
of these models is that thermonuclear burning begins 
as a subsonic, turbulent deflagration and then makes 
a transition to a supersonic, shock-driven detonation 
(Khokhlov 1991ab, Yamaoka et al. 1992, Woosley & 
Weaver 1994). These models arc generally known as 
delayed detonation models. In this class of models the 
amount of nickel produced is a function of the den- 
sity at which the transition is made from deflagration 
to detonation, the central density, metallicity and the 
chemical structure of the inital WD (Hoflich et al. 
1995, Hoflich & Khokhlov 1996, Hoflich, Wheeler & 
Thielemann 1998, HWT98 hereafter). The radioac- 
tive decay of the variable nickel mass is the dominant 
factor which gives a range in maximum brightness and 
^^Ni is the most important factor which governs the 
light curve shape. The models with less nickel are 
not only dimmer, but are cooler and have lower opac- 
ity, giving them redder, more steeply declining light 
curves. A given amount of nickel can be produced by 
different combinations of the model parameters and, 
from the models, we expect a spread around the mean 
maximum brightness-decline relation of w 0.3 to 
0.5™. A similar spread of k, ±0.4" {a « 0.18" in 
B ) in the local, observed relation between maximum 
brightness and decline (Hamuy et al. 1996) may also 
suggest that there are additional parameters needed 
in the empirical relations. Note that new observations 
and recalibrations of the old observations indicate a 
somewhat tighter relation (« 0.12™, Jha et al. 1999; 
Kirshner et al., 1999). From the current status of 
theoretical models, this very narrow spread cannot be 
understood, but it cannot be ruled out either. Both 
reanalysis of new light curves and theoretical investi- 
gation of the coupling between the "free" parameters 
may help to answer this puzzle. Note that this also 
implies that observable spreads in relations, e.g. be- 
tween rise time and decline time, cannot be expected 
to change along the path given by the change of a 
particular free model parameter. 

There are some hints that SN la have undergone 
evolutionary effects. Branch et al. (1996) have shown 
that the mean peak brightness is dimmer in ellipti- 
cals than in spiral galaxies. Wang, Hoflich & Wheeler 
(1997) found that the peak brightness in the outer 
region of spirals is similar to those found in ellip- 
ticals, but that in the central region both intrinsi- 
cally brighter and dimmer SNe la occur. This implies 
that the underlying progenitor populations are more 
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inhomogcncous in the inner parts of spiral galaxies 
which contain both young and old progenitors. Sys- 
tematic effects that must be taken into account in the 
use of SNe la to determine cosmological parameters 
include technical problems, changes of the environ- 
ment with time, changes in the statistical properties 
of the SNe la, and changes in the physical properties 
of SNe la. 

In a previous paper (HWT98), the effect of the 
change of the initial metallicity and the C/0 ratio 
on light curves and spectra was studied by chang- 
ing both the C/0 ratio and the metallicity indepen- 
dently. It was found that a reduction in the C/0 
ratio mainly effects the energetics of the explosion. 
The ^^Ni production are reduced and the Si-rich lay- 
ers are more confined in velocity space for smaller 
C/0. Moreover, the brightness to decline ratio can 
change causing a shift in the zero point of the max- 
imum brightness decline relation of « 0.3™. Such a 
shift can be detected by a corresponding change in 
the rise time by about 3 days. The influence of the 
metallicity is especially important with respect to the 
^^Fe production in the outer layers which govern the 
spectra in the UV and blue; however, the influence on 
the resulting rest frame visual and blue light curves 
is found to be small. A separation of the effects of Z 
and the C/0 ratio in the WD is justified because the 
metallicity is inherited from the ISM when it is formed 
and the C and O are produced during the progenitor 
evolution. To first order, the C/0 ratio depends on 
the zero age main sequence mass. Consequently, a 
systematic shift in the maximum brightness/decline 
relation can be produced by a change with redshift of 
the typical progenitor mass. This effect is expected as 
the mean life time becomes shorter if we are looking 
back closer to the star forming period. 

To second order, however, one must take into ac- 
count that the metallicity also influences the stellar 
evolution during central helium burning. The possi- 
ble consequences for evolutionary effects in SNIa and 
the influence of ^^C(q!, 7)-rate as has been pointed 
out recently by various groups (Nomoto et al. 1997; 
Dominguez et al. 1999; Hoflich et al. 1999; Umeda 
et al. 1999a). Thus, besides changing the C/0 ratio 
as a consequence of the change in the typical progen- 
itor mass, the effects of metallicity discussed above 
may also apply to individual SNe la of a given main 
sequence mass. 

In this paper we study the latter effect: the in- 
fluence of the metallicity on an individual supernova 



with respect to the nucleosynthesis and light curves. 
We investigate the size of the shift of the zero point in 
the brightness/decline relation and the consequences 
for the use of SNIe la for cosmology. We address the 
issue how to distinguish the effect of a change in the 
"typical" progenitor from the effect of the metallicity 
on the chemical structure of an individual progenitor. 

2. Brief Description of the Numerical Meth- 
ods 

2.1. Stellar Evolution 

The stellar evolution has been calculated using the 
code of Nomoto's group up to the end of the he- 
lium burning. These calculations have already been 
published and discussed in detail by Umeda et al. 
(1999a). Subsequently, the evolution of the C/0 core 
is calculated by accreting H/He rich material at a 
given constant accretion rate on the core by solv- 
ing the standard equations for stellar evolution us- 
ing a Henye scheme. Nomoto"s equation of state is 
used (Nomoto et al. 1982). Crystallization is ne- 
glected. For the energy transport, conduction (Itoh 
et al. 1983), convection in the mixing length the- 
ory, and radiation are taken into account. Radiative 
opacities for free-free and bound-free transitions are 
treated in Kramer's approximation and by free elec- 
trons. A nuclear network of 35 species up to '^'^Mg is 
taken into account. Here, we do not include convec- 
tive mixing at the onset of ignition of the deflagration 
front at the center. Depending on the description, the 
entire WD may mixed be mixed (Nomoto, Sugimoto 
& Neo 1976, Hoflich et al. 1998). Although this may 
change the central carbon concentration at the center, 
the energetics of the explosions will not change. 

2.2. Hydrodynamics 

The explosions are calculated using a one-dimensional 
radiation-hydro code, including nuclear networks (Hoflich 
& Khokhlov 1996 and references therein). This code 
solves the hydrodynamical equations explicitly by the 
piecewise parabolic method (Collela & Woodward 
1984) and includes the solution of the frequency aver- 
aged radiation transport implicitly via moment equa- 
tions, expansion opacities (see below), and a detailed 
equation of state. Nuclear burning is taken into ac- 
count using a network which has been tested in many 
explosive environments (see Thielemann et al. 1996 
and references therein). 
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2.3. Light Curves 

Based on the explosion models, the subsequent ex- 
pansion and bolometric as well as broad band light 
curves are calculated using a scheme recently de- 
veloped, tested and widely applied to SN la (e.g. 
HWT98 and references therein). The code used in 
this phase is similar to that described above, but nu- 
clear burning is neglected and 7 ray transport is in- 
cluded via a Monte Carlo scheme. In order to al- 
low for a more consistent treatment of scattering, 
wc solve both the (two lowest) timc-dcipendent, fre- 
quency averaged radiation moment equations for the 
radiation energy and the radiation flux, and a total 
energy equation. At each time step, we then use 
T(r) to determine the Eddington factors and mean 
opacities by solving the frequency-dependent radia- 
tion transport equation in the comoving frame and 
integrate to obtain the frequency-averaged quantities. 
About one thousand frequencies (in one hundred fre- 
quency groups) and about five hundred depth points 
are used. The averaged opacities are calculated un- 
der the assumption of local thermodynamical equi- 
librium. Both the monochromatic and mean opaci- 
ties are calculated using the Sobolev approximation 
(Sobolev 1957). The scattering, photon redistribu- 
tion and thermalization terms used in the light curve 
opacity calculation are calibrated with NLTE calcula- 
tions using the formalism of the equivalent-two-level 
approach (Hoflich 1995). 

3. Results 

The initial mass of the C/0 WD is given by the 
results of stellar evolution. Its mass dc;pends on Mms 
of the progenitor and Z. At the time of the explosion, 
the WD masses are close to the Chandrasekhar limit. 
The WD has grown by accretion of H/He and subse- 
quent burning. In the accreted layers, the C/0- ratio 
is close to 1. 

Here, we study the size of the metallicity effect us- 
ing as an example a 7 Mq with Population I (Z=0.02) 
and II (Z=0.004) compositions (Fig. 1). We assume 
that the relative abundance within the metals is solar. 
In reality, the [O /Fe] changes with the metallicity and 
we may separate between Zcno and Zpe- The for- 
mer governs the CNO cycle including the production 
of neon which determines the proton to nucleon ratio 
and, consequently, the abundances of the iron group 
elements produced during the thermonuclear explo- 
sion. This will effect the spectra at maximum light 



(see HWT98) but hardly the size of the C/0 ratio 
which is the main subject of this letter. It is the iron 
abundance that governs the opacity and changes the 
core size. The metallicity Z of heavy elements, e.g. 
Fe, mainly effects the convection during stellar He- 
lium burning and, consequently, the size of the C/0 
core and the central C/0 ratio. We note that the ex- 
act size of the effect and its sign depends on the mass 
of the progenitor on the zero age main sequence, and 
Z. Even for a given mass, the changes are not mono- 
tonic, but may change sign from Pop I to Pop II to 
Pop HI (Umeda et al. 1999a, Dominguez et al. 1999). 
In addition, the variation depends sensitively on the 
assumed physics such as the ^^C(a, 7)^^0-rate (e.g. 
Straniero et al. 1997). Our example can serve as a 
guide to estimate the size of this effect, but will not 
give all possible variations. 

At the time of the explosion, the central density 
has been fixed at pc = 2AE9g/cm~^. The ratio a of 
the deflagration velocity to the sound speed has been 
set to 0.02 and the transition density ptr at which the 
deflagration is assumed to turn into a detonation is 
is 2.4 lO''^ cm~^). The parameters are close to those 
which reproduce both the spectra and light curves 
reasonably well (Nomoto et al. 1984; Hoflich 1995; 
Hoflich & Khokhlov 1996). 

The total abundances of the most important ele- 
ments are given for the model in Table 1. Fig. 2 
gives the density and velocity versus mass for the 
models with the Pop I and Pop II metallicity, re- 
spectively. Fig. 3 gives the composition proflles of 
the major elements for the same models. The over- 
all density and velocity structures are insensitive to 
changes in the underlying model. Qualitatively, the 
flnal burning products can be understood by the re- 
lation between hydrodynamic and the nuclear time 
scales. Since most of the energy is released by the 
explosive burning of carbon and oxygen, the energy 
production per gram depends on the initial chemical 
composition of the WD, i.e. the C/0 ratio. The nu- 
clear time scales are determined by the peak temper- 
ature during burning, which depends on the energy 
release per volume because the energy density is ra- 
diation dominated. Since the energy release per gram 
is fixed; the peak temperature is given by the compo- 
sition and the local density. Thus, the latter variables 
are the dominant factors that determine the final com- 
position of a zone. After accretion on the initial core, 
the total mass ratio of Mc/Mq is 0.75 and 0.61 for 
the models with Pop I and II, respectively. The total 



4 



Table 1: Total abundances of the delayed detonation model for high and low metallicity 



Population 


He 


C 





Mg 


Si 


S 


Ca 


Fe 


I 
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Fig. 1. — Chemical profile of the C/0 core of a star 
with the main sequence mass of 7Mq after central 
Helium burning. 
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Fig. 2. — Density and velocity as a function of mass 
for the delayed detonation model originating from 
progenitors with PopI and Pop II composition. 
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Fig. 3. — Same as Fig. 2, but abundances as a func- 
tion of the final expansion velocity. Both the initial 
^^Ni and the final Fe profiles are shown. 



nuclear energy release depends on the total C/0 ratio, 
but the structure of the WD is unaffected by a change 
in metaUicity (HWT98). Therefore, the final kinetic 
energy is reduced from 1.41a;10^^erg for Z=0.02 to 
1.37xl0^^erg for Z=0.004. The nuclear energy pro- 
duction during burning decreases with the C/0 ratio 
at the burning front and, consequently, the transi- 
tion density is reached later in time. The resulting 
larger pre-expansion of the outer layers reduces the 
'^^Ni production from 0.56 to 0.51 M© for the Pop I 
and II model, respectively, and the transition regions 
between complete and incomplete C and O burning 
changes by about 500 km s^. 

As discussed in HWT98, a change of the metallic- 
ity changes also the abundances of trace elements such 
as ^''Fe (Fig. 4). The reason is that the metallicity 
mainly affects the initial CNO abundances of a star. 
These are converted during the pre-explosion stel- 
lar evolution to "'^'^N in H-burning and via ^^N(q;,7) 
i8p(^+) i8o(c,, j) 22]s^e to nuclei with N=Z+2 in He- 
burning. The result is that increasing metallicity 
yields a smaller proton to nucleon ratio, Yg, through- 
out the pre-explosive WD (Thielemann et al. 1997, 
HWT98). Higher metallicity and smaller Yg lead to 
the production of more neutron-rich Fe group nuclei 
and less ^^Ni. For lower metallicity and thus higher 
Ye, some additional ^^Ni is produced at the expense of 
^^Fe and ^^Ni (Thielemann, Nomoto & Yokoi, 1986). 
We note that these layers with v> 12,000 km 
dominate the spectra around maximum light and lines 
of the iron group, but have little influence on the light 
curves. 

Light curves provide a valuable tool to probe the 
underlying explosion models, namely the absolute 
amount of ^^Ni and its distribution. Broad-band light 
curves are shown in Fig. 5. As expected from the 
last section, the main effect on the light curves is 
caused by the change in the explosion energy and 
the ^^Ni production. The change in the maximum 
brightness remains small (My (PopI) — My (PopI I) = 
—0.03™) and the rise times are different by about 1 
day {tB[[PopI) = 16.4d vs. tB{PopII) = 17. 6d). 
The smaller expansion due to the smaller Ekin causes 
a reduced geometrical dilution of the matter and a re- 
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Fig. 4. — Same as Fig. 2, but abundances of different 
isotopes as a function of the expansion velocity for 
layers with partial burning. 
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Fig. 5. — Comparison of light curves in B and V of 
the delayed detonation models with metallicities cor- 
responding to Pop. I and Pop. II. 



duction of the expansion work at a given time. Con- 
sequently, for Pop II, the rise time to maximum light 
is slower by about 1 days and the maximum bright- 
ness is only slightly smaller despite the 10 % reduced 
production of ^^Ni because, in the latter case, more 
of the stored energy goes into the radiation rather 
than kinetic energy, and energy can be stored for a 
longer time. The most significant effect is a steeper 
decline ratio and a reduced ^®Ni production for the 
Pop Il-model. This translates into a systematic off- 
set of w 0.1™ in the maximum brightness decline ra- 
tio somewhat depending on the slope derive from the 
data (Phillips 1992, Hamuy et al. 1996). Using ei- 
ther the stretching method or the LCS-method gives 
similar offsets. 

4. Conclusions 

The influence of the metallicity at the main se- 
quence on the chemical structure of the exploding 
white dwarf, the nucleosynthesis during the explosion, 
and the light curves of individual Type la supernovae 
have been studied using the example of delayed deto- 
nation models Q Model parameters are used for which 
the theoretical optical and infrared light curves and 
of the spectral evolution resembles those of a normal 
bright SNe la. Note that we did not tune the param- 
eter that the specific example resembles a particular 
supernova or a specific normalization such as a stretch 
factor s=l. 

A change of the metallicity influences the central 
helium burning and, consequently, the size of the 
C/0 core which becomes a C/0 WD and its C/0 
ratio. The WD may grow to the Chandrasekhar 
mass and explode as a Type la supernovae. Con- 
sequently, the metaUicity changes the C/0 structure 
of the exploding white dwarf. As C and O are the fuel 



^Addendum: In this paper, we consider the dependence of the 
brightness-decline rate relation and changes in the rise times on 
Xc for given model parameters. We assume that the C/O ra- 
tio has little effect on the nuclear burning. In contrast, Umeda 
et al. (1999b) considered the influence of Xc on the transi- 
tion density and, thus, on the production of radioactive Ni. 
Fully consistent with previous papers (Hofiich 1995, Hoflich, 
Khokhlov & Wheeler 1995), they find that the Ni produc- 
tion declines with the transition density between deflagration 
and detonation. The absolute brightness declines accordingly 
(HK96). Umeda et al. (1999) do not consider the effect of 
Xc on the brightness-decline rate relation or the rise times. 
Instead, for an assumed relation between transition density 
and C/O ratio, they discuss the implications on the typical 
brightness as a function of red-shift based on their study for 
progenitor evolution. 
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for the; thermonuclear explosion, the metallicity indi- 
rectly changes the energetics of the explosion. In ad- 
dition, the metallicity alters the isotopic composition 
of the outer layers of the ejecta that have undergone 
explosive O burning. Especially important is the in- 
crease of the ^^Fe production with metallicity which 
alters the spectra near maximum light (HWT98). 

As the C/0 ratio of the WD is decreased, the 
explosion energy and the ^^Ni production are re- 
duced. Changing the initial metallicity from Pop I 
to Pop II changes the rise times by about 1 day, and 
causes a small decrease in luminosity at maximum 
light, a faster post-maximum decline and a larger ra- 
tio between maximum light and ^^Ni tail by about 
0.1™. This effect is equivalent to a change from a 
one-parameter maximum brightness/decline relation 
which is widely used to determine the cosmological 
parameters Om and Oa- In our example, the calibra- 
tion of the AM(rfMi5) relation is changed by about 
0.1™. From previous studies of the stellar evolution 
for different masses and metallicites ( Umeda et al. 
1999a, Dominguez et al. 1999), the influence of Z on 
the C/0 ratio is comparable to our example, but the 
sign may change depending on details of the nuclear 
cross sections and the treatment of convection, and 
the mass and metallicity range. In our example, the 
systematic change in the progenitor metallicity from 
Pop I to Pop II translates into a shift of the confidence 
regions in the Qm — ^^A plane based on high-z SNel a 
(Riess et al. 1998; Perlmutter et al. 1999) by about 
0.15 and 0.7 along the small and large axes of the er- 
ror ellipsoids. However, a change of typical progenitor 
from Pop I to II is unlikely for the redshifts currently 
used to determine and Cl\ because we see a large 
variation of metallicities in our galaxy (Edvardsson et 
al., 1993ab). 

Even strong variations in the metallicity will hardly 
produce selection effects with respect to the counting 
rates for SNe la at large red-shifts. Variations of the 
expected size are critical if, in future, SNe la may be 
used to measure large scale scalar-fields because the 
metallicity may show large local variations during the 
early phases when individual SNe explosions govern 
the metallicity. 

A change of the C/0 ratio reveals itself mainly by 
the change in the ratio between rise time and decline. 
This change may be caused either by a change of Z or, 
alternatively, by a change of the main sequence mass 
of the progenitor. When using statistical methods to 
determine cosmological parameters, both a change in 



the distribution of progenitor masses and the metal- 
licity of the sample may cause systematic evolution- 
ary effects. Fortunately, it does not matter whether a 
change of C/0 is due to a change in the metallicity or 
the typical mass of the progenitors. In both cases, we 
find that AM » Atrise where AM is the offset in the 
maximum brightness/decline relation in magnitudes 
and At is measured in days (compare HWT98). 

The effect of progenitor mass and metallicity can 
be untangled by simultaneous analysis of both spectra 
and light-curves. In principle, this allows a method 
to get a handle on the progenitor mass. With the 
current level of modeling such analysis is restricted 
to differential comparisons between individual super- 
novae. 

After submission of this paper, a comparison of rise 

times between local and distant supernovae has been 
submitted (Riess et al. 1999). They find a difference 
in the rise times between the local and the distant 
sample of 2.5 ± 0.4 days if all observed supernovae are 
normalized to the stretch factor s of 1. This result 
is based on a preliminary analysis of the high red- 
shift data by Goldhaber (1998). It may also be noted 
that both the local and the distant sample show an 
spread of ±lday in the rise time/decline relation with- 
out overlap which is hard to understand. However, 
if confirmed and in light of our analysis, this would 
imply a systematic offset w 0.25™ in the brightness 
decline relation. A change of this order would not al- 
ter the basic conclusion of the high z searches that we 
live in a a low O universe but, then, a low fl universe 
with Hm = 0.2 and Oa = may be consistent with 
SN data as well. 

The following trend for the theoretical models is 
worth noting: For realistic cores, the both the mean 
Mc/Mo and the Mc/Mq in the central regions of the 
WD tend to be smaller than the canonical value of 1 
used in all calculations prior to 1998 (e.g. Nomoto 
et al. 1984, Woosley & Weaver 1994, Hoflich & 
Khokhlov 1996). Consequently, as a general trend, 
the rise times are about 1-4 days slower compared to 
models published before 1998. 
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